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Abstract. 

A novel realization of the classical SU(2) algebra is introduced for the Dirac 
relativistic hydrogen atom defining a set of operators that, besides, allow the 
factorization of the problem. An extra phase is needed as a new variable in 
order to define the algebra. We take advantage of the operators to solve the 
Dirac equation using algebraic methods. To acomplish this, a similar path to 
the one used in the angular momentum case is used; hence, the radial eigenfun- 
tions calculated comprise non unitary representations of the algebra. One of the 
interesting properties of such non unitary representations is that they are not 
labeled by integer nor by half-integer numbers as happens in the usual angular 
momentum representation. 
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I. Introduction 

The unitary representations of groups and algebras are of a great interest 
in physics. We can mention the ubiquitous example of angular momentum or 
the SU(2) algebra. In this instance, as it is well known, the representations 
are labeled by two real numbers j and m which may take only integer or half- 
integer values; the representations depend only on two parameters, 9 and <f> say, 
which are defined in the compact sets [0, it] and [0, 2tt] respectively. 1 However, 
there are many other physically interesting groups or algebras which are not 
unitary. In such a case, the representations are not restricted to parameters 
defined in compact sets nor its generators are necessarily Hermitian, but they 
nevertheless can play an important role in physics. The Lorentz group being 
a very important example of a physically relevant non compact group whose 
algebra is not necessarily unitary. 2 

It is the purpose of this article to introduce a realization of the cyclic SU(2) 
algebra in terms of non Hermitian operators and then to use these operators to 
factorize and solve the relativistic Dirac hydrogen atom. The solution is obtained 
using algebraic methods using the basic operators of the system, following a route 
parallel to the unitary compact case, but introducing an extra variable (found 
below to play the role of a phase) which is required in our approach. 3 

The use of algebraic techniques has been common for these sort of problems. 
For example, in the relativistic hydrogen atom a la Dirac one approach has 
been the use of shift operators; 4 whereas in the non relativistic case a succesful 
approach requires the use of ladder operators and the factorization method. 5-6 
Our approach is more akin to the introduction of ladder operators than to the 
shift operator method customarily used for this problem. 4 

The paper is organized as follows. In section II we introduce the equations 
of the problem and define our notation. In section III we construct the basic 
operators spanning the SU(2) algebra as a useful tool for the problem. In section 
IV we define the inner product needed to investigate the properties of our basic 
operators. In section V some non unitary representations of the symmetry al- 
gebra are constructed as the radial eigenfunctions of the problem. We find that 
they are non compact representations that play an equivalent role to the spheri- 
cal harmonics in the unitary case. In order to write such representations we find 
convenient to define a family of polynomials that are found to be the associated 
generalized Laguerre polynomials. In section VI, using algebraic methods and 
the operators defined in section II, we find the energy spectrum of the hydrogen 
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atom. In section VII we give our conclusions. In the Appendix, we list and plot 
the explicit expressions of the first 6 mentioned polynomials; they are explicitly 
related to the generalized associated Laguerre polynomials previously used by 
Davis for expressing the radial eigenfunctions of the hydrogen atom. 7 

II. The Dirac hydrogen atom 

The symmetry algebra for the bound states of the non-relativistic hydro- 
gen atom, even in the classical case, is well known to be the SO(4), 8 ' 9 but in 
the relativistic case the situation is different as it is well known. The algebra 
associated with the symmetry of the radial part of the problem can be regarded 
as SU(2) as occurs with the angular momentum. The beautiful thing is that 
we can then proceed to solve the relativistic hydrogen atom following a method 
that essentially parallels the calculation of eigenfunctions and eigenvalues of the 
angular momentum at the only expense of introducing an extra phase. 

Let us begin with the Dirac Hamiltonian of the hydrogen atom 

Ze 2 

H D =a-p + f3m , (1) 

r 

where m is the mass of an electron and a and (3 are the standard Dirac matrices 10 

-(:;)■ o-(s !)■ (2) 

where the l's and O's stand, respectively, for 2 x 2 unit and zero matrices and the 
a is the standard vector composed by the three Pauli matrices a = (a x , a y ,a z ). 
Since the Hamiltonian (1) is invariant under rotations, the solutions of the prob- 
lem can be written in the form 



V(r,M) = - 
r 



/ F(r)y jm (6, 



(3) 

\iG{r)y jm {e,<j>) " 

Where 3^ m and y'j m are spinor spherical harmonics of opposed parity. Parity 
is a good quantum number in the problem because the Coulomb potential is 
invariant under reflections; parity goes as (— I) 1 and, according to the triangle's 
rule of additon of momenta, the orbital angular momentum / is given by Z = j± \. 
But, instead of working directly with parity, we prefer to introduce the quantum 
number e defined by 



(4) 
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Thus I = j + | in all cases; we also define I' = j — |. Accordingly, the spherical 
spinor y^ m depends on / whereas the spherical spinor 3^ m , which has an opposite 
parity, depends on /'. Writing the solutions in the form (3) completely solves the 
angular part of the problem. 

Let us now address the radial part of the problem; we are interested in its 
bound states, then the quantity k = \Jm 2 — E 2 is positive definite; furthermore, 
let us define 



<-^, r jSeW + i), -^f, (5) 

then, we can write the differential equations for the radial part of the problem, 
in terms of the dimensionless variable p = kr, as 

v -4z + -)g(p)=(-" + -)f(p), (6) 

and 



dp p J V P 



these equations are to be regarded as the initial formulation of our problem. 

The first thing we want to do is to show that Eqs. (6) and (7) can be 
rewritten using a set of three operators whose commutation relations define a 
SU(2) algebra. To this end, let us first introduce the new variable x through 
the relation (but please notice that this change is not required for any of the 
conclusions that follow, see the second remark below after Eq. (15)) 

P = e x , (8) 

so x is defined in the open interval (—00,00) and redefine the radial functions 
F(p) and G(p) : introduced in equations (6) and (7), in the form 

F(p(x)) = V^+E ty_(x) + Mx)} , (9) 



G(p(x)) = Vm-E [if>-(x) - i)+(x)} . (10) 

In terms of the new functions ^>+(x) and ^_(x), we thus arrive to the 
following set of equations for our problem 
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and 



*- {i)= h^w +T ) Mx) - (12) 

This first-order system can be uncoupled multiplying by the left the first equation 
(Eq. (11)) times the operators that appear between square brackets in the second 
equation and, viceversa, by multiplying the second equation (Eq. (12)) times the 
operators that appear (between square brackets) in the first one. This procedure 
gives us the second order system 



d 

— -e x + 
ax 



(E 



E 2 



+ 2pe x - e 2x - 



dx A 



ip + {x) 



(13) 



and 



ip-{x) 



(14) 



where we have defined 



fj, 



+ 1. (15) 

ym 2 — E 2 

At this point there are several remarks that need to be done. First, in the 
next section the seemingly odd term —1/4 in Eqs. (13) and (14) is shown to 
be necessary to construct the algebra. Second, as we said before, the change 
of variable from p to x is not really necessary for any of the calculations in the 
article, however, we prefer to work in the x rather than in the p variable be- 
cause this choice simplifies the appearance of some of the equations and, mainly, 
because it makes the inner product introduced in section III (in Eq. 33) look 
more familiar; but to make contact with the usual description we sometimes, at 
our convenience, revert to the variable p. As a third remark, notice that we can 
regard equations (13) and (14) as two eigenvalue equations where the common 
eigenvalue u> is given by 



u J = rf-C 2 - 1 -=j(j + l)-e; (16) 

as it is obvious, we do not need to calculate u> because it follows directly from 
the radial symmetry of the problem and from the intensity of the interaction 
which is needed to set the scale. The fourth remark we want to do is that, as 
the minimum value of j is 1/2, then u > for at least Z = 1, 2, • • • up to 118; 
for a discussion of the significance of this number see Ref. 4, page 236. 

III. An operator algebra for the Dirac hydrogen atom 

As the main purpose of this article is the construction of non unitary rep- 
resentations of the SU(2) algebra for the Dirac hydrogen atom, 3 let us first 
introduce the operator 

^3 = ~i§- e (17) 

depending exclusively on a new variable £, which is essentially an extra phase as 
must be clear in what follows, and the two operators 

which depend both on £ and on the transformed 'radial' variable x. These three 
operators satisfy the following algebraic relations 

[fi 3 ,0±] = ±Q±, (19) 

and 

[0+,0_]=20 3 . (20) 
We can alternatively define the two operators Oi and 2 as 

1 = ^(0 + + 0_), 2 = l(0+-0_), (21) 
in terms of which the algebraic properties (19) and (20) read 

[Oi, Qj] = ieijkQk, i,j,k = 1,2,3. (22) 

It is now clear that the conmutation relations (19) and (20) — or just (22) - 
correspond to an SU(2) algebra. 1 ' 3 ' 8 To complete the discussion, we also need to 
introduce the Casimir operator of the algebra; let us consider the operator 
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o 2 = vt\ + o 2 + n 2 , 



(23) 



we can easily show that (23) is indeed a Casimir for the algebra (19) and (20) 
(or (22)) 

[fi 2 ,fij = 0. for % = 1,2,3. (24) 

Here, as in the angular momentum case, we can regard O 2 as the square of 
Q = fi]i + 2 j + 3 k. To obtain the explicit expression of O 2 , it is better to 
calculate first the product 0_0 + : 

8 2 8 8 8 2 1 

«-«+^-^- 2 '4 +i l + 5H (25) 

and obtain the Casimir operator from the relationship O 2 = + 03(^3 — 1). 

We remark that we do not have a linear term in d/dx in equation (22) because 
we choose the last constant in equation (18) precisely as 1/2. We then easily 
conclude that the Casimir operator is given explicitly by 

^ = ^-e--2^|-i. ,26) 

Although we are not restricted to a compact set of parameters anymore due 
to the presence of the variable x, in analogy with the conventions used for the 
spherical harmonics Y" z m (6>, (/>), we label the simultaneous eigenfunctions of O 2 
and 3 as V£(x,a); where the numbers u are the eigenvalues of O 2 and the 
numbers fi are the eigenvalues of f^. Therefore, we write 



W(x,0 = /*^(a;,0 (27) 

and 

n 2 v»(x,o=oov»(x,0; (28) 

we thus have 

V£{x i S) = e i '*V£{x). (29) 

Where, again, we have used the angular momentum analogy to write the x 
functions as V£(x), resembling the usual notation for the Legendre polinomials 
P[ n (9). In this equation it becomes clear the role of £ as just an extra phase. 
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Notice that Eqs. (26) and (28) serve to establish the connection between 
the introduced operators of the SU(2) algebra and the radial part of the hy- 
drogen atom problem, since applying the Casimir operator O 2 to V£(x,£) or to 
V^ -1 (#,£), reproduces Eqs. (13) and (14). Furthermore, comparison with these 
same equations tell us that 

lMs) = P£(z), (30) 

^(x)=Vr\x). (31) 

The operators Q± play the role of ladder operators for the problem; they move 
along the set of eigenfunctions changing the eigenvalue fx to the eigenvalue fj, ± 1, 
in a completely analogous way to the the case of the angular momentum algebra, 

n ± v^(x,o^v^ ±1 ( x ,o, (32) 

the required proportionality constants are evaluated in section V (Eq. (52)). 

IV. The inner product and related properties of the operators 

To establish the properties of the O-operators and to construct the rep- 
resentations of the SU(2) algebra they span, we need a properly defined inner 
product. Here we cannot longer rely on the angular momentum analogy, since 
two of the parameters of the algebra we purport to construct are defined over the 
non-compact interval (—00,00), making it completely different from the angular 
momentum. This situation implies that not all the generators of the algebra 
are to be Hermitian (indeed, O3 is found to be Hermitian but fii and O2 are 
are found to be anti-Hermitian) and, as a consequence, that the corresponding 
group could not be unitary. 

Let us consider the operation 

(0,VO = y ^ J_ <f(t,x)if,(t,x)dx. (33) 

It is easy to show that Eq. (33) defines an inner product, since, as it is not 
difficult to prove, it satisfies the three basic properties: 

i ) If 1>(x) = 0, then (ip,tp) = 0. 
ii ) If ij){x) ^ 0, then (V>, > 0. 

in ) If c is any complex number, then (-0, c0) = c(-i/>, 4>) and (cip, 4>) = c*(ip, 4>). 
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To study the behavior of the generators of the algebra in terms of the inner 
product (33), we consider a certain function ip(x, £) associated to a certain fixed 
value hq. We define the elements of the associated Hilbert space H^ with 
functions of the form 

ip(x, = e l ^ 0+m - n) ^{x), (34) 

where m and n are integer numbers and F(x) is a well behaved function depend- 
ing only on x. As we exhibit in section V, this is actually the general form of 
the functions inhabiting the Hilbert space of our problem. 

With the help of the inner product (33), we can establish the properties of 
the O a (a = 1, 2, 3) operators. Let us consider first fi 3 , in this case the important 
product is 

/oo 
T*{x)T{x)dx; (35) 
-oo 

this last equation, due to the presence of the Kronecker delta, can be written as 



f 2n dp f°° r b 

/ — / i^- e -^ 0+m '- n '^ e- l ^ 0+rn - n) ^*(x)T(x)dx, (36) 
Jo 2tt J_ 00 I d£ 



and this is precisely (O3?//, ip). We have in this way proved that O3 is an Her- 
mitian operator: 

fij = 3 . (37) 
To study the operators 0±, let us first consider fi + and evaluate the product 



(v>',ft+vo = r — r e-^ o+m '- n,) ^*(x)dx 

JO 27T J^ 

e%i ~ eX ~ + I) ei(f " 0+m ~ n) ^ :F ( x ) dx - 



(38) 



To analyse this integral, it is simpler to split it in two parts. Let us consider first 
the e^(— <9/<9£ + 1/2) part; its contribution to the inner product (38) is 

^io + m-n+^j 6 m '-n',m-n+i J J 7 * {x)J 7 (x) dx, (39) 
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this expression can be written as 



2tt 



_a i 

d£ + 2 



J(fi +m' -n' 



(40) 



e <(Ai0+m - n)€ ^(x) dx. 



We now take care of the contribution of the term e^(<9/<9£ — e x ). After a partial 
integration this contribution becomes 



W-n'-l,m-n 



<9x 



it 



+ e 1 :F(x) 



JF(x) dx. 



(41) 



Taking together the two previous results, it is easy to see that the operator 
complies with = — fi_. A similar calculation establish the analogous property 
fit = — fi + . Therefore, we have established that 



0± = — O =f . (42) 

We have proved that not all the operators are Hermitian thence arriving 
to the anticipated results: First, we showed that the operator fi 3 is indeed 
Hermitian, as a consequence we expect the range of the only parameter it depends 
on, £, to be compact; this is certainly the case since £ G [0,27r], a compact set. 
The other two operators fii and Vt 2 depend upon x, a variable defined over the 
non compact set (— oo, oo), but here such operators are anti- Hermitian 



Qt=_Q a , a = 1,2 (43) 
as it is easy to show from Eqs. (21) and (42). 

V. Representations of the SU(2) algebra (radial eigenfunctions of the 
problem). 

With the inner product defined in the previous section, we are in the position 
of introducing a complete orthogonal basis of simultaneous eigenfunctions for 
O 2 and 3 which accordingly must carry a representation of the algebra — and, 
besides, they solve our radial eigenvalue problem (Eqs. (6) and (7)). We have 
decided to choose the conmuting operators in similar fashion to the standard 
SU(2) case. Let us define then 



V^(x^) = \ufi>, 
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(44) 



where the kets \u> n > are both assumed orthogonal and normalized respect the 
inner product (33) 

< (jj'n'\uH >= <L,u;'<W- ( 45 ) 

Now, as Oi and Q2 are not Hermitian, the Casimir operator Q 2 defined in 
(23) is not necessarily positive definite, we can nevertheless introduce a positive 
definite operator as 

O f O = -O 2 - fi 2 , + = 2Qj - O 2 ; (46) 
the positivity of this operator allows us to show that 

2/j 2 > iv, (47) 

that is, is bounded by below. As a consequence, there must exist a minimum 
value for let us say X = m j n . We also know the that the ket \u X > is 
annihilated by fi_ or, equivalently, that Q + Q-\uj X >= 0; so 

w _A 2 -A = or oj = X(X-1). (48) 

Given equation (47), let us introduce a slight change in notation, writing A 
instead of uj in the eigenfunctions V%(x) of Eq. (29), that is, V%(x) is to be 
replaced by V%(x). Furthermore, as A has to be positive and since uj = r 2 — ( 2 — 
1/4, we find that the minimum -value is 

Mmin = A = s +^ ( 49 ) 

where we have defined 

s^+sJr]-e = \J(j + l/2) 2 -e. (50) 

Notice that s is a real quantity for Z = 1, 2, 3, ... , 137. 

However, the most important conclusion we can draw from Eqs. (49) and 
(50), is that A no longer has to be an integer or half-integer number, as necessarily 
happens in the unitary compact angular momentum case. We may see that this 
result is a direct consequence of both fix and not being Hermitian operators. 
So the operators introduced for the SU(2) algebra associated with the Dirac 
hydrogen atom lead naturally to non unitary representations labeled by real 
numbers, or, at least, not necessarily integer nor half-integer ones. 
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For constructing the functions comprising the representations, let us intro- 
duce the constants as follows (compare with Eq. (32)) 

Q±\ufi >= C^\u ^ ± 1 >; (51) 

these constants can be explicitly evaluated from < u fi\Q + Q,- \u> (i >= C~C^_ X 
and from (C^_i)* = if we further assume these constants to be real, we 

easily get 

C± = ± V / M^±1)-A(A-1), (52) 

a result with a slightly different form than in the analogous angular momentum 

case. 4 ' 8 

The ground state of the hydrogen atom, can be obtained from the equation 
fi_|AA >= for the positive set of eigenvalues. Such equation becomes 



(53) 



whose solution is 

V${x) = d A e sx exp(-e x ) = d x e (A " 1/2)x exp(-e x ), (54) 

where 

2 (A-l/2) 

lh s 7m=T) <56) 

is a normalization constant and T(y) stands for the Euler-gamma function. Since 
A is the lowest eigenvalue, in this instance we should have 

= V£(x), (56) 

and 

ip_(x) = Q. (57) 

In terms of the radial variable p, we can see from Eqs. (9) and (10) that the 
ground state solutions are 
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and 



G( P ) = - 



(m-E) 



-p 



(59) 



2m(A - 1/2) 



normalized in the sense 




(60) 



We can also see that in the case of negative eigenvalues the solutions behave 
as ~ p s e p , giving divergent behavior as p — > oo. This behavior makes the negative 
energy solutions not square integrable; therefore, we have to discard them if we 
want to describe physically realizable states. 2 ' 3 ' 10 

The excited states can be obtained, as we do in the appendix, by applying 
succesively Q + to the ground state | A A >, as it is customarily done for the 
spherical harmonics. The result involves the functions V%(x) introduced by 
Eqs. (29) and (34). The V%(x) are polynomials multiplied by the weight factor 
W(p) = p x - x l 2 e-P. This weight factor assures that the behavior of the big 
and the small components of the spinor are regular both at the origin as well 
as p — > oo. As illustration of the solutions discussed here, we quote in the 
Appendix the first few cases of V x +P (x), for p = 0, . . . , 5, there we also plot 
these polynomial part of the functions (i. e. we plot V%(x) without the weight 
factors) for the first 5 excited energy levels. Notice that the polynomial part in 
the eigenfunction V^ +P (x) is always of order p. 

The matrix representations of the Q a , a = 1,2, 3, for X = s + 1/2 + p, p = 
0, 1, 2, 3, . . ., are constructed from equations (27), (28), (51) and (52). In contrast 
with the standard SU(2) case, here the representations are non Hermitian and 
infinite dimensional for each A, excepting for O3. The matrix elements of the 
operator O3, including the negative eigenvalue series, are given by 



where p = ±(A + p), p = 0, 1, 2, . . .; therefore, its trace vanishes and the deter- 
minant of any element of the corresponding group with the form exp(zfi3)£, is 
always 1. 

For the other two operators the only non-vanishing matrix elements are 



< up\Qs\ujp' >= p5 j 



MM'' 



(61) 



< up\Vti\ujp ± 1 >= =F- a/M/" ± 1) - MS ~ 1)> 



(62) 
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< w//|f2 2 |u>/x ± 1 >= --y/nbj, ± 1) - A(A - 1), (63) 

This means again that the trace of both Oi and fi 2 vanish. Notice also that the 
determinant of a group element generated by or fi 2 is 1 only for a purely 
imaginary parameter. 

VI. The energy spectrum. 

We can now evaluate the bound energy spectrum for the problem. As we 
mentioned in section V, the bound state energy spectrum comprises only the 
positive series eigenvalues. 2 ' 10 Let us first express the energy in terms of the 
eigenvalue \x from Eq. (15) 

-1/2 

(64) 

then, for the case of positive eigenvalues — the only ones with physically ap- 
propriate eigenfunctions — we have that \x = A + p, where p is an non negative 
integer p = 0, 1, 2, . . ., or, equivalently, that \i — 1/2 = s + p. This gives precisely 
the energy spectrum of the relativistic hydrogen atom. To rewrite our result in 
a more familiar form, we need only to define the principal quantum number n 
and the auxiliary quantity ej as follows 

n = j + ^+p, (65) 

e i = n - s - p = 3 + \ - s ; ( 66 ) 

we then finally conclude that /z — 1/2 = s + p = n — 6j, which gives precisely the 
well-known energy spectrum. 10 

VII. Concluding remarks. 

In summary, we have constructed an SU(2) algebra for the hydrogen atom 
in the Dirac formulation, introducing the Hermitian operator and the anti- 
Hermitian operators fii and fi 2 . The result of all of this, is that the represen- 
tations are labelled by numbers A wich are neither integers nor half-integers as 
in the case of the more familiar unitary representations. Nevertheless, the al- 
gebra introduced predicts precisely the energy eigenvalues and eigenfunctions of 
the Dirac hydrogen atom. One of the most noteworthy features of the repre- 
sentations reported here is the mixing of a spinorial angular momentum char- 
acter, implying an equally spaced spectrum, with the energy requirements of 



E 



m 



1 + 



(jjl- 1/2)2 
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the problem — requiring a differently spaced spectrum; the interplay of these two 
spectral requirements is basically reflected in the fact that the eigenvalues associ- 
ated with Eqs. (13) and (14) follows from both the generic radial symmetry and 
the specific features of the interaction in the Dirac equation (1). From Eq. (15) it 
also follows that in the limit of vanishing interaction, i. e. £ — > 0, our representa- 
tions collapse and, in this special case, \i = 1 always. Such behavior is precisely 
as expected because there is no longer any restriction over the eigenvalues and 
thus the spectrum becomes continuous, corresponding to a free Dirac particle. 
The operator algebra we introduced allows an essentially algebraic solutiond of 
Dirac hydrogen atom which may have various applications. 3 ' 4 ' 6 

It is to be noted also the possible connections that our formulation may have 
with systems with hidden supersymmetric properties, 13-16 as we will discuss in a 
forthcoming article. The energy spectrum of the problem has some peculiarities 
which also appear in the spectrum of a Dirac oscillator; 14-15 ' 17 in particular the 
equally spaced energy solutios for if)+(x) and ip-(x) resemble, respectively, the 
behavior of big and the small components of the aforementioned system. As a 
result of this resemblance, we are studying the hidden supersymmetric properties 
and the superconformal algebra, in the sense of Refs. 15 and 18, associated with 
this problem. The similitude of Eqs. (13) and (14) with the corresponding ones 
for a Morse oscillator should be also noticed. 6 

It is worth pinpointing that we are forced to introduce the new variable £ 
in order to define the algebra; in terms of the solutions of the Dirac equation, 
£ just plays the role of a phase. To exemplify, when we perform a "rotation" 
using 3 , the big component changes from F(x) oc [ip_(x) + tp+(x)] to F(x) oc 
[ e *(M-i)^^_( x ) + e^ifj+ix)] = e v*[ e -i€iMaO+iMaO]. The phase e*^ does not 
play any observable role, but the term e~ 1 ^ changes the relative phase between 
the ip+(x) and the ip-(x) components of the eigenfunction and, in consequence, 
changes the radial function F(p) itself although the energy spectrum is still 
invariant under such transformation; this is a consequence of the fact that \oj \x > 
and Qs\uj n > both correspond to the same eigenvalue fx. In a way, this resembles 
what happens when there are superselection rules in a system. 19-21 

To finalize, let us comment that it is not widely known that the radial eigen- 
functions of the Dirac hydrogen atom can be expressed in terms of generalized 
associated Laguerre polynomials, as was realized by Davis a long time ago. 7 ' 11 
These polynomials, which are a generalization to non integral indices of the usual 
associated Laguerre polynomials, are defined as 7 
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c ^ x) = T ^rt P +i) lFl{ - p;a + 1;x) (67) 

where T(x) is again the Euler gamma function, p is a positive integer, and the 
iFi(—p, a+1; x) stands for the confluent hypergeometric function — having one of 
their arguments negative, the hypergeometric function reduces to a polynomial. 12 
The polynomial representation of the radial eigenfunctions introduced here is 
related to that used by Davis in Eqs. (A8) and (A9) of the Appendix. 
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Appendix 

The purpose of this Appendix is to illustrate the behavior of the first few 
polynomials associated to the radial eigenfunction of the hydrogen atom in the 
base | A p > and to relate this description to the old but little known results of 
Davis. 7 ' 11 We explicitly calculate the first 6 functions of the positive eigenvalue 
series, such functions are always of the form the weight factor W(p) times a 
polynomial; the weight factor is W(p) = p s e~ p . The polynomials are plotted 
in Figure 1. The radial eigenfunctions are essentially generated from the basic 
relationship Q + \uj p >= C+\u p + 1 >. In the p variable, the first equation in 
the series can be written as 

e+li { p ifp~ p ~ l ik + e ^ (p) = c M +el(M+1) ^A A+1 (p) m 

which is just the first term in the whole ascending series V^ +P (p) = O^P^(p) 
used to recursively calculate (A2-A7). Note that the polynomial part of the 
function V^ +P (x) is always of the form J2 P =0 Ci(X)p\ where Cj(A) is also an 
order (p — i) polynomial in A. 

Starting with p = A, the first few functions in the positive series are then 
given by 



16 



V X M =W(p), (A2) 
r^ +1 (p)= ] fl(\-p)W(p), (A3) 

p ^ M= vmrw) l2p2 - (2X+1)i2p - x)] ' {M) 

+A(2A 2 + 3A + 1)], A(5) 

V*+ 4 (p) =4/1 - W{p) = [p 4 - 2 (3 + 2A)p 3 + 3 (3 + 5A + 2A V 

A 1 V 3 ^/A(A + 1/2)(A + 1)(A + 3/2) U V " V " 

- (3 + 11A + 12A 2 + 4A 3 ) p + (A 3 + 3A 2 + ^A + ^)A], (A6) 

X v/15A(A + 1/2)(A + 1)(A + 3/2) (A + 2) L H 2 JI 

- (10A 2 + 51A + 54) p s + (10A 3 + 66A 2 + ^A + ^) p 2 

— ■ — ■ 

, rl4 319 2 257 x 33, x5 15 x4 65 x3 

- 5A 4 + 37A 3 + —A 2 + —A + — ) p + X + —A 4 + —A 3 

4 4 2 2 4 

105, 2 27 , 
+— A 2 + -A]. (A7) 

The relationship to these polynomials to those used by Davis can be seen 
from Eqs. (13) and (14), putting p = expx and introducing the function v(p) 
according to ip + = p s exp(— p) v(p), and finally using C{p) = v(p/2), to get 



d r / _ . «, -■ d£> 

p df +[{2s + 1) - p] Tp + 



(s 2 + C 2 -rf) 



1/2) 



£ = 0, 



048) 



where /it = s + l/2+p. This equation can be regarded as a generalization to non 
integer index of the usual associated Laguerre differential equation. The equation 
corresponding to ip- can be obtained in analogous fashion. This equation reduces 
to the one used by Davis if and only if s 2 = t 2 — ( 2 , a result which just recovers 
the definition in Eq. (50) . Now we can give the explicit relationship between our 
representation of the radial eigenfunctions to that used by Davis as 



iMp) =v 



_ T) s+l/2+p 



(p)=p s exp(-p)£ 2s (2p), 



s+1/2 

■/: 
/s 



(AS) 
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save for normalization factors (unimportant for the point at hand) where s = 
A — 1/2, and where the generalized associated Laguerre polynomials used by 
Davis are defined in Eq. (67). 
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Figure Caption 



Figure 1 

We show the graph of the first polynomials V^ +p (p)/W(p), for p = 1,2, 3, 4, 5, 
Z = 1 and j = 1/2. Notice the similarity of the behavior of all polynomials near 
the origin. 
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